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ABSTRACT 



We present a comparison of radio observations of the archetypal powerful 
radio galaxy Cygnus A and 2-D numerical hydrodynamical simulations. We 
characterize some global trends in the observed radio properties and compare 
them with the properties of a simulated radio source. The observational results 
are the following. The width of the observed surface brightness distribution per- 
pendicular to the source axis can be well characterized by a Gaussian over most 
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of the length of the source. The ratio of the Gaussian FWHM to the second 
moment is fairly constant along the source with an average value of about 2.5 
indicating that they give roughly consistent measurements of the source width. 
The average observed surface brightness, estimated pressure, and estimated min- 
imum energy B field decrease with distance from the hot spots. We find evidence 
for significant structure in the estimated cross-sectional slices of emissivity. The 
numerical results are the following. Jets propagating in a constant density at- 
mosphere will decelerate with time. Thus, the estimated dynamical age of the 
source will be greater than the actual age of the source. For a source similar to 
Cygnus A the difference will be about a factor of 2. The second moment gives 
an accurate representation of the "true" width of the simulated source. The 
Gaussian FWHM tends to be about 40% larger than the true width and can 
be systematically in error if the surface brightness exhibits multiple peaks. We 
suggest that the ratio of the Gaussian FWHM to the second moment may be a 
diagnostic of the emissivity profile in the lobes. The simulations can qualitatively 
reproduce the overall observed morphology and the behavior of the cross-sections 
in surface brightness, the decline in surface brightness with distance from the hot 
spots, and the width of the lobes. This suggests that the 2-D simulations give a 
reasonable representation of the properties of Cygnus A. 

Subject headings: galaxies: jets - galaxies: intergalactic medium - hydrodynam- 
ics - radio continuum: galaxies 



1. Introduction 

The primary physical processes in large, powerful FRH radio galaxies are still poorly 
understood. There is a consensus that the sources are powered by two oppositely directed, 
highly collimated outflows from an AGN - resulting in shocks and radio emission (Scheuer 
1974; Blandford & Rees 1974). It is clear that the youngest populations of relativistic 
electrons are located near the radio hotspots at the extremities of the source (e.g., Leahy, 
Muxlow, & Stephens 1989; Carilli et al. 1991; Alexander & Leahy 1987). Beyond this, 
differences of opinions begin to emerge. 

It is important to understand the physical processes that shape the observed radio 
emission of a source. If these could be identified and understood, it might be possible to 
use the observed source properties to determine the pressure, density, and temperature of 
the ambient gas, and the beam power of the source. These could then be used to probe the 
evolution of source environments with redshift, and evolution of AGN properties, such as jet 
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power. These, in turn, would suggest the circumstances that lead to AGN activity, and the 
role of environment in both triggering (perhaps via mergers or close encounters) and then 
constraining (through the gaseous environment) the radio properties of a source. 

Much effort has been spent on analysis of multi-frequency radio observations of radio 
galaxies (e.g., as summarized in Leahy 1991), and independently on numerical simulations 
of radio galaxies (e.g., Norman et al. 1982; Lind et al. 1989; Clarke, Norman & Burns 1989; 
Hardee & Norman 1990; Falle 1991; Clarke & Burns 1991; Cioffi & Blondin 1992; Loken 
et al. 1992; Mioduszewski, Hughes & Duncan 1997; Marti et al. 1997; Komissarov & Falle 
1997, 1998; Rosen et al. 1999; Carvalho & O'Dea 2002a,b). 

Here we attempt to bridge these two fields by directly comparing radio images of the 
archetypal radio galaxy Cygnus A with the results of 2-D hydrodynamical simulations of a 
source which is expected to be similar to Cygnus A. Our goal is to explore the behavior of 
global observed properties and relate these to intrinsic properties of the underlying flows. 



2. Radio Observations & Data Analysis 

Two observations of Cygnus A are studied in detail. The 151 MHz radio image from 
Leahy, Muxlow, & Stephens (1989) was kindly provided to us by Paddy Leahy, and the 1345 
MHz image from Carilli et al. (1991), was kindly provided to us by Chris Carilli. The 151 
MHz image, obtained with MERLIN, has an angular resolution of 3 arcsec down and mean 
and rms noise levels of 0.1 and 0.3 Jy/beam. The 1345 MHz image, obtained with the VLA, 
has an angular resolution of 1.36 arcsec down to mean and rms noise levels of 0.01 and 0.01 
Jy/beam. 



2.1. Surface Brightness Profiles and Gaussian Widths 

The FITS images of the source at 151 and 1345 MHz were rotated in an identical 
manner so the hot spots lie on a horizontal line. The hot spot-hot spot axis was taken to be 
a symmetry axis of the source. Using AIPS and auxiliary programs, vertical slices of surface 
brightness perpendicular to the hot spot symmetry axis were taken at one pixel intervals of 
the FITS image, as illustrated in Figure 1. Each surface brightness profile was then studied 
in detail. Distances to each slice are measured in arc seconds relative to the hot spot on each 
side of the source, increasing as we move in toward the source center. 

Each slice profile was fit with a Gaussian curve to determine the FWHM of the best 
fitting Gaussian. In these fits, we included the radio data above a threshold of the mean plus 
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three times the rms noise level of each image; data below this level were removed. These 
cuts were taken at 1 Jy/beam for the 151 MHz image and 0.04 Jy/beam for the 1345 MHz 
image. 

The radio surface brightness profiles and best fit Gaussian curves are shown at integer 
intervals of 3 arc seconds (the FWHM of the 151 MHz image), in Figures 2, 3, 4, 5 for the 
151 and 1345 MHz data, and for each side of the source. 

Generally, the Gaussian fits provide a good fit to the profiles; occasionally the edges of 
the profiles drop more steeply than the wings of the Gaussian fit. 

The full width at half maximum of these Gaussian fits are shown as a function of slice 
distance from the hot spot at 151 and 1345 MHz in Figure 6. The errors on the Gaussian 
widths are much smaller than the FWHM of the radio observations, so the one sigma error 
bar on the FWHM of the Gaussian was taken to be half of the FWHM of the observing 
beam. Thus, an uncertainty of 1.5 arc seconds for the 151 MHz data, and 0.7 arc seconds 
for the 1345 MHz data was adopted as the uncertainty of the Gaussian widths. 

The ratio of the best fit Gaussian widths at 151 and 1345 MHz are shown in Figure 
7. The fact that this ratio is unity over most of the source is reassuring, and suggests that 
to first order the 151 and 1345 MHz data define similar source shape and structure. The 
average value of the ratio for the left and right hand sides of the source, respectively, are 
W^i5i/W^i345 = 0.94 ± 0.01 and 0.90 ± 0.01 when all the data are included; when only data 
with hot spot distances between 6" and 48" are included these average value are 0.92 ± 0.01 
and 1.01 ± 0.01 for the left and right sides of the source, respectively. 



2.2. The First and Second Moments 

Two other measures of source structure are the first and second moments of the cross- 
sectional surface brightness profiles described above. These were obtained using the standard 
expressions x = '^{xi 5'j)/^(S'j), and \/^{Si [xi — x]"^)/ ^{Si), where Si is the surface 
brightness at the point Xi, shown as the surface brightness and x-axis in Figures 2 to 5. 
These are plotted as a function of the slice distance from the hot spot in Figures 8 and 9. 
The uncertainty of the first moment is taken to be half the beam width, as above. 

The first moment is a measure of the location of the surface brightness weighted central 
axis of the source. In the 151 and 1345 MHz data, we see a varying displacement (or 
wandering) of the surface brightness weighted source axis which is small compared with the 
Gaussian FWHM of the bridge at similar bridge locations. The first moment can be affected 
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by the jet that is apparent in the radio grey scale images, and it is interesting that the 
deviations of the first moment from zero (expected for a perfectly cylindrically symmetric 
source) are in opposite sense for the right and left-hand side of the source, which could 
be influenced in part by any wobbling of the outflow axis of the jets. A detailed analysis 
indicates, however, that there is not an exact correspondence between the deviations of the 
first moment on each side of the source. 

The second moment is another measure of the bridge width, and is shown in Figure 9. 
The bridge width as measured by the best fit Gaussian and the second moment are compared 
in Figure 10, which shows the ratio on the best fit Gaussian FWHM to the second moment. 
For the 151 MHz data, this ratio is 2.59 ± 0.07 for the left side of the source and 2.49 ± 0.07 
for the right side of the source, while for the 1345 MHz data, this ratio is 2.45 ± 0.03 and 
2.53 ± 0.03 for the left and right hand sides respectively. The fact that the ratio is similar 
over most of the source suggests that either method may be used to estimate the bridge 
width. The physical interpretation of the numerical value of the ratio is discussed in §3 
(numerical simulations). Note that the ratio Wc/W2 discussed in §3.4, obtained from the 
numerical simulations, is half of the value stated here since W2 is twice the bridge radius 
indicated by the second moment. 



2.3. The Average Surface Brightness Along the Bridge 

The average surface brightness of each cross-sectional surface brightness profile is shown 
in Figure 11. This is obtained by summing up the total area under a particular surface 
brightness profile (see Figures 2 to 5) and dividing by the FWHM of the best fit Gaussian 
curve, Wg'- S = '^{SiAx)/WG, where Si is the surface brightness for the ith point. Ax is 
the interval along x around the point i. Thus, the units of the average surface brightness are 
Jy/beam. Numerical simulations (see §3) indicate that Wq is a good estimate of the bridge 
diameter, suggesting that the surface brightness values obtained and shown in Figure 11 do 
not need to be re-normalized. 

The uncertainty of the average surface brightness is obtained by adding the rms noise 
level of the image in quadrature with the beam width: {6S/S)^ = (rms/S)^ + {Ob/^Wc)^, 
where Ob is the FWHM of the observing beam; recall that the uncertainty on the Gaussian 
width is taken to be half the beam size. 

The contrast between the hot spots and the bridge is much higher for the 1345 MHz 
data than it is for the 151 MHz data, and the values of S differ substantially for the data 
sets. But, once we move beyond the hot spots into the bridge region, the profile of the 
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average surface brightness along the bridge axis is similar for the two frequencies, which can 
be seen in Figure 11. 



2.4. The Average Pressure and Minimum Energy Magnetic Field Along the 

Bridge 

The average bridge pressure for each slice is plotted as a function of distance from the 
hot spot on each side of the source in Figure 12. This is a "minimum energy" pressure 
(Burbidge 1956) in the sense that it is obtained using the equation P = (S/Wg)^^'^ (see 
Figures 2 to 5). Thus, the pressure has units of (Jy/arc sec)^/^. To convert this to physical 
units of erg/cm^ and thus obtain the normalization factor to apply to the pressure presented, 
we note that the pressure is simply related to the minimum energy magnetic field (e.g. Wan 
et al. 2000): P = (1.336^^'^ + 6^)-B^/(247r), where B is the minimum energy magnetic 
field strength and b parameterizes the offset from minimum energy conditions: the true field 
strength is bB. As described below, the minimum energy magnetic field strength estimated 
using the 151 MHz data has a normalization factor of 30/iG, so the normalization factor for 
the pressure estimated using the 151 MHz data is about 1.2 x 10^^^ erg cm~'^(1.336^^'^ + 6^), 
or about 2.8 x 10~^^ erg cm~^ for 6 = 1, or about 1.3 x 10"^'^ erg cm~^ for b = 0.25, as 
suggested by Carilli et al. (1991), and confirmed by Wellman, Daly, & Wan (1997a,b) for 
the bridge region of Cygnus A. At 1345 MHz the normalization factor for the pressure is 
]^g-io g^g cm~^(1.336~^-^ + 6^), or about 2.4 x 10~^° erg cm~^ for 6 = 1, or about 1.1 x 
10^^ erg cm~^ for b = 0.25. 

It is clear from Figure 12 that the minimum energy magnetic fields computed for the 
bridge region using the 151 MHz and 1345 MHz data are similar in magnitude and structure, 
while the pressures computed for the hotspots are considerably different - with the 1345 MHz 
data indicating a much higher pressure than the 151 MHz data. 

The average minimum energy magnetic field strength is each slice is obtained by taking 
the square root of the minimum energy pressure, B = {S/WgY^'^ . The normalization factor 
is obtained using the standard expression for the minimum energy magnetic field strength, 
assuming zero energy in relativistic protons, a filling factor of unity, a constant spectral 
index of unity, frequency cut offs of 10 MHz and 100 GHz, and a cosmology with a value 
of Hubble's constant of 70 km/s/Mpc, a normalized mean mass density today of Qm = 0.3, 
and a cosmological constant of Q\ = 0.7. These parameters yielded a normalization for 
the minimum energy magnetic field strength for the 151 MHz data of about 30/iG, and a 
normalization factor of about 88/iG for the 1345 MHz data. The minimum energy magnetic 
field strengths are shown in Figure 13. Like the minimum energy pressure, the minimum 
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energy magnetic field strengths across the bridge computed using the 151 MHz and 1345 
MHz data are quite similar in magnitude and structure, while the hotspot magnetic field 
indicated by the 1345 MHz data is much larger than that indicated by the 151 MHz data. 
One of the reasons that the estimated pressures in the radio hot spots at 151 MHz and 1345 
MHz differ is that the magnetic field strengths and hence pressures have been estimated 
assuming a constant spectral index of unity. While this may be a reasonable estimate in the 
bridge region of the source, it clearly is not a good approximation in the hot spot region. 

Note that the total pressure P in relativistic electrons and fields estimated assuming 
that this pressure is proportional to the square of the minimum energy magnetic field implies 
that P oc (e)^^^, where e is the volume emissivity of the radio emission (see §2.5). This follows 
from the fact that the radio surface brigthness S is obtained by integrating the emissivity 
along the source path length, so S* ~ eWc- Thus, P = {S /WqY^'^ , which follows directly 
from the minimum energy argument, implies that P oc e^/^. This is a good approximation 
whenever the true field strength is proportional to the minimum energy magnetic field with 
a proportionality constant that is independent of position in the source. 



2.5. The Radio Emissivity as a Function of Position 

The emissivity, or energy emitted per unit volume per unit time at one particular fre- 
quency can be obtained from the surface brightness profiles shown in Figures 2 to 5. For a 
particular surface brightness profile, it is assumed that the emission is cylindrically symmet- 
ric, and that the diameter of the cylinder is equal to the total width of that profile (only data 
above the noise level of l*mean + 3 rms are included). Then, the emissivity as a function 
of radius from the center of the slice can be determined. This is done in the following way. 
The data point that is furthest from the center of the slice with a surface brightness above 
the noise level is used to determine the edge of the source. We move in an amount Ax to 
the position that is centered on i?i, and compute the emissivity Ci at the point i?i by taking 
the surface brightness at this point divided by the length of the chord through the circle 
that that point, as illustrated in Figure 14. We then have the emissivity ei at all points 
that are at this distance, i?i, from the center of the cylinder for that slice. Stepping in, 
the next surface brightness data point can be used to determine the emissivity €2 that is a 
distance R2 from the center of the circle; this is done by accounting for the fact that the 
total surface brightness 5*2 will have contributions from ei and €2, illustrated in Figure 14. 
This process is continued to obtain the emissivity as each distance i?„ from the center of the 
axis of symmetry of the source, and is valid to the extent that the slice possesses cylindrical 
symmetry. 



The normalization for the volume emissivity e^ = dE/{dV dt dv) can be obtained by 
noting that the emission coefficient jy = e,^/(47r) for a volume element that emits isotropically, 
and the specific intensity I^, = ju^y, where Ay is the length of the ray through the source 
that is shown in Figure 14. The observed specific intensity I^^ = Ii,{l + z)~^, and has units 
of Jy/beam. Using these expressions it is easy to show that the normalization factor for the 
emissivities shown in Figures 15, 16, 17, and 18 are 1.8 x 10~^^ erg s~^ cm~^ Hz~^ for the 
151 MHz data, and 9.0 x lO"^'^ erg s"^ cm^^ Hz~^ 

These emissivities show significant structure, indicating that in any given slice there can 
be regions with very high and/or very low emissivities (or that the assumption of cylindrical 
symmetry was not a valid assumption for that slice). A relatively large change in emissivity 
is needed to produce a relatively small surface brightness feature. Some of the changes in 
emissivity are suggestive of a channel or jet-like feature, which are suggested by either an 
emission feature, or a cavity in the emissivity near the center of the slice. The variations 
in emissivity could be due to variations from cylindrical symmetry, hydromagnetic waves, 
turbulence, or other effects. 



3. Cygnus A versus Numerical Simulations 

Here we present the results of 2-D numerical hydro simulations and compare them with 
the observational results for Cygnus A. The advantage is that in the simulations we know 
the intrinsic properties of the "source" . The disadvantage is that the simulations may not 
include all the relevant physics which occurs in real radio galaxies. 



3.1. The Numerical Simulations 

We have carried out axisymmetrical hydrodynamical numerical simulations of light jets 
propagating in a constant density atmosphere. The details of the simulations are discussed 
by Carvalho & O'Dea (2002a,b). We have chosen not to include magnetic fields in the 
simulations. If the magnetic field is dominant, MHD simulations produce structures (e.g., 
the "nose cone", Clarke et al 1986; Lind et al. 1989) which do not appear to correspond 
to the observed radio structures in Cygnus A. If the magnetic field is not dominant, the 
overall structure and dynamics of the source in the MHD simulation are similar to the case 
where the magnetic field is absent . In addition, in a particle pressure dominated source the 
magnetic field distribution closely follows the distribution of gas pressure and density (Lind 
et al 1989). Numerical simulations by Komissarov (1989) show that if the magnetic field is 
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initially weak, it will come into equilibrium with the gas pressure near the cocoon boundary 
because the gas pressure drives the sideways expansion of the cocoon, and the gas is not able 
to cross the magnetic field lines. 

Numerical simulations studies (e.g., Norman et al. 1982; Carvalho & O'Dea 2002a,b; 
Krause 2003) indicate that light jets will infiate a cocoon whose width is inversely propor- 
tional to the density ratio [r]) between the jet and the ambient gas. The over expanded 
bridge in Cyg A implies a very low value for t] (< 10^^). According to Alexander and Pooley 
(1996) an order-of-magnitude estimate of the jet density gives r/ ~ 4 x 10~^. We have run a 
number of simulations where we used 77 = 2 x 10~^ and a jet internal Mach number M = 10 
(Run DOl), r]= Ix 10-^ M = 11 (Run D02) and 77 = 1 x 10-^ M = 10 (Run DOS). One 
of the reasons for doing so is to ensure that small changes on the jet parameters do not have 
a significant effect on the final results. Since the overall linear size of Cyg A is ~ 120 kpc 
we used a simulation grid with physical dimension 66 kpc in the direction of the source axis 
and a jet radius 1 kpc. We expect that, as the bow shock in front of the jet head reaches 
the end of the grid, the expanded cocoon would give a fair representation of the source lobe. 
We have run the simulations with several different grid widths (equivalent to 16.5, 33 and 
49.5 kpc) in order to study the effect on the cocoon of the bow shock leaving the grid. In 
the largest grid (Run DOl), the entire source is contained on the grid, while in the smaller 
grids, the bow shock eventually leaves the grid laterally. We find that the details of the 
cocoon properties differ in these simulations, though the overall properties of the cocoons 
are similar. This gives us confidence in the results of simulations carried out in the narrower 
grids which are computationally much cheaper. 

The final stage of the simulation of Run DOl is shown as a density contour map in Fig. 
19. After comparing this with Runs D02 and D03 we concluded that no significant change 
is observed in the source physical properties although their general appearance may look 
different. Therefore, in what follows we show the results for Run DOl only. In Fig. 20a we 
show the contours of the bow shock, cocoon and jet. Fig. 20 also shows the average pressure 
inside the cocoon (b) the lateral expansion speed of the cocoon (measured at the contact 
discontinuity) (c) and the lateral expansion speed of the bow shock (measured at the outer 
shock front) (d) along the source axis. 

The average pressure was obtained in the following way. For each slice, with width equal 
to the grid cell and having cylindrical symmetry, the total thermal energy was calculated by 
summing up the product of the energy density per cell times the volume element of the cell, 
and then dividing by the total volume of the slice to obtain the total internal energy density, 
0.5pf^„g or 3/2nkT. This was then multiplied by (F — 1) to obtain the average pressure in 
that cell. The simulations was carried out for a monotonic, ideal, non-relativistic fluid, thus 
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ratio of specific heats is F = 5/3. 

The shock structure seen in the numerical simulation is shown in Figure 20a. As usual, 
three shock regions can be identified: the outer shock front, the contact discontinuity, and 
the inner shock, which includes the entire bridge region. 

In the simulation, the location of the jet can also be identified. The outer shock front is 
also called the bow shock, and is indicated by the outermost line in Figure 20a. The contact 
discontinuity separates the shocked ambient gas from the shocked jet fluid and is indicated 
by the middle line in Figure 20a. The radio emitting fluid that would produce the observed 
radio emission is expected to lie within the contact discontinuity. Thus, the region interior 
to the middle line in Figure 20a is what is referred to as the "cocoon" throughout this paper. 
The shocked ambient gas lies between the bow shock and the cocoon and is not discussed 
further in this paper. The jet is identified in the simulations by means of trace particles 
which are injected together with the jet gas. 

Slices are taken across the cocoon, and the average pressure within the cocoon is com- 
puted for each slice as described above. This average pressure is plotted as a function of 
distance from the hot spot in Figure 20b, and is normalized to be unity at the location of 
the hot spot. This measure of the pressure includes only the thermal pressure of the "gas" 
within the cocoon of the numerical simulation and does not include kinetic pressure. In the 
numerical simulation, the average pressure in slices across the cocoon undergo a series of 
oscillations as we move along the symmetry axis of the bridge. 

The lateral expansion of the cocoon (i.e., as defined by the contact discontinuity) and 
that of the bow shock are shown in Figures 20c and d, and are normalized to be unity at the 
location of the hot spot. Like the bridge width, and the average pressure in a slice across 
the bridge (i.e., the cocoon), the lateral expansion speed of the cocoon undergoes a series 
of oscillations rather than varying smoothly along the contact discontinuity. This variation 
in lateral expansion speed is reflected in differences in the bridge width. For example, in 
the simulation the bridge width is seen to have a bulge at distances from the hot spot of 
about 10 to 30 kpc, and this corresponds to a region where there is an increase in the lateral 
expansion speed of the cocoon relative to nearby points. The lateral expansion speed of the 
bow shock varies much more smoothly and monotonically with distance from the hot spot. 
And, appears to generally follow Ub oc Z^^/^, which is predicted theoretically for the regions 
in which the expansion speed exceeds the ambient sound speed, after which the expansion 
speed is expected to level off (e.g., Begelman & Ciofii 1989; Daly 1990, 1994; Wellman et al. 
1997a; Reynolds, Heinz, & Begelman 2002; Carvalho & O'Dea 2002a). 



-11- 



3.2. Dynamical and Synchrotron Age 

The ratio between the source advance speed and the initial jet speed is shown in Fig. 
21 as a function of source size. For comparison, we also show the ambient sound speed 
(dashed line). We see how the source decelerates rapidly with its Mach number relative to 
the ambient sound speed decreasing by a factor of ~ 22 from ~ 70 to ~ 3.2. This is due to 
the progressive increase of the jet head radius. The observed speed fluctuations are also due 
to oscillations of the head radius. These oscillations of the head radius are probably driven 
by vortex shedding (e.g., Lind et al 1989; Norman, Winkler, & Smarr 1982). We note that 
changes in jet momentum flux may also drive changes in the advance speed. Near the end 
of the simulation, an increase is observed in the advance speed which experience indicates is 
a transient phenomenon. 

Two effects likely to occur in nature that are not included in the simulation may also 
come into play. In nature, the density of the ambient gas may be slowly decreasing with 
distance from the central engine rather than remaining constant, and this would cause the 
decline in the head velocity to be less significant. On the other hand, in nature the outflow 
axis of the jet may "wobble" rather than pointing in a flxed direction (as in the simulation), 
which would go in the opposite sense, and at any given time would cause the velocity of the 
head to be lower (e.g. Scheuer 1982; Cox, Gull, & Scheuer 1991). This second effect could 
also modify the bridge width as a function of distance from the hot spot. 

The dynamical age of Cygnus A determined using the ram pressure of the hotspot is 
about an order of magnitude larger than the estimated synchrotron age (Carilli et al. 1991). 
It has been suggested by some authors that this discrepancy can be alleviated if one relaxes 
the equipartition hypothesis when estimating the magnetic fleld. According to Carilli et al. 
the magnetic fleld must be 1/3 of the equipartition value. Similar results have been found by 
Perley and Taylor (1991) for 3C295, and by Wellman, Daly and Wan (1997a,b) for a large 
sample of sources. The deceleration of the source head observed in the simulations can, at 
least in part, be responsible for the discrepancy. 

Thus, the age calculated using the speed at the source maximum length is larger than 
the real age of the source. In our simulation the difference is not too large because at the 
end of the simulation the head speed has momentarily increased. However, if we calculate 
the age by taking the head speed when the source was, for instance, ~ 51 kpc in size, it 
would be larger than the real age by a factor ~ 2. Thus, this effect can account for only 
a part of the discrepancy between the ram-pressure balance and synchrotron ages found by 
CariUi et al. (1991). 
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3.3. Surface Brightness and Radio Emissivity 

We have calculated the volume emissivity of the source assuming that the energy den- 
sity of the relativistic electron component is proportional to the hydrodynamical (thermal) 
pressure of the gas inside the cocoon and an equipartition magnetic field. This is easy to 
do since the radio emissivity is proportional to P^/^, as described at the end of §2.4. We 
will also consider separately contributions proportional to the thermal pressure, the kinetic 
pressure, and the total pressure (see Figure 25). We assume that the radio emission comes 
from the jet itself and the cocoon (Leahy, 1991) where the pressure is higher than ambient 
pressure. This region (corresponding to the region interior to the middle line in Figure 20a) 
is filled with jet material which has been shocked in the head and expands laterally. Besides 
any non-relativistic gas that may be present in this region, it is assumed that the region also 
contains the basic ingredients necessary to generate synchrotron emission, that is, a mixture 
of magnetic field and relativistic electrons most probably accelerated by Fermi process at the 
strong shock region in the jet head. We integrated the radio emissivity along the line of sight 
(assuming the source lies in the plane of the sky) to obtain the surface brightness, which 
allows us to draw synthetic radio images of the source. We have also taken into account 
radiation loss of the relativistic electron population in a very simple way. We divided the 
cocoon into slices and assigned an age to each slice equal to the time that has passed since 
it has been laid down by the advancing jet head. The modified spectrum is then calculated 
using the Kardashev-Pacholczyk model (Kardashev 1962; Pacholczyk 1970) supposing that 
the magnetic field is the equipartition field. 

In this manner, we were able to predict the radio structure that would be observed 
from the source indicated by the numerical simulation. In Fig. 22 we show the predicted 
151 MHz radio image of the source at five evolutionary stages; the synthetic data have been 
smoothed by a Gaussian with FWHM of 3 arcsec to match the 151 MHz data, as shown 
by the black circle in panel e of the figure. The age of the source corresponding to each 
panel is: (a) 1.35 Myr, (b) 4.53 Myr, (c) 6.17 Myr, (d) 8.89 Myr and (e) 9.15 Myr. One 
can clearly see the characteristic edge-brightened shape of powerful FR-II sources and an 
aspect ratio resembling that of Cygnus A. When the source is young, radiative energy losses 
are negligible and we can see the radio bridge extending all the way from the head to the 
central galaxy (panels b and c). As the source ages, the central parts becomes less luminous 
and it assumes the classical double lobe shape (panel e). The hotspots are evident in all the 
simulated images. 

Cross-sectional slices of the surface brightness distribution for the source as it appears 
in panel e, when it is about 9 Myr old, are shown in Fig. 23. The second panel in the first 
row corresponds to the position of the hotspot while the bottom ones show the parts of the 
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source near the central object. We have normahzed the surface brightness to unity at the 
hot spot and the scale in each row is the same as that in Fig. 4 at 151 MHz. The simulated 
surface brightness is qualitatively similar to that observed in Cygnus A (Figs. 2, 3, 4, 5). Of 
course, the simulation is axisymmetric and thus so are the cross-sectional surface brightness 
slices, whereas in nature this is not always the case. 

The Gaussian fit is shown as a dotted line in the figure. The central spikes that appear 
in some of the panels are caused by the jet contribution to the surface brightness. The spikes 
are more evident in those regions of the jet where it has been compressed (Fig. 19). The 
Gaussian fits provide a reasonable approximation to the actual surface brightness slices over 
most of the source except along row 3 where the wings of the Gaussian fit clearly extend 
well beyond the boundary of the source. 

Fig. 24 shows the predicted average surface brightness of a slice as a function of distance 
from the hotspot in the numerical simulation; the average surface brightness is obtained as 
described in section 2.3, but using the true width of the simulated source. 

Here we investigate how the different forms of energy (thermal pressure, kinetic pressure, 
and total pressure) could contribute to the radio emission. Again, the pressure in relativistic 
electrons and magnetic fields is assumed to be proportional to the gas pressure, this being 
determined by the local thermal pressure, kinetic pressure, or total pressure given by the 
numerical simulation. The thermal pressure is the thermodynamic pressure associated with 
the gas microscopic motion (discussed in §3.1). The kinetic pressure is equal to (F — 1) 
times the kinetic energy density, where the kinetic energy density 0.5pf^ depends on the 
large-scale, macroscopic or bulk velocity v of the gas. The total pressure is the sum of 
the thermal and kinetic pressures. We also use the approximation that the emissivity is 
proportional to P'^/^ (described in §2.4). Considering the kinetic pressure allows for the fact 
that macroscopic kinetic energy in the form of turbulent motion could be transferred to 
the relativistic component through the Fermi process or through hydromagnetic instabilities 
(e.g. De Young 1980; Eilek & Henriksen 1984; Eilek & Shore 1989). For this figure, the radio 
emissivity is calculated supposing that the local pressure of the relativistic fiuid, including 
the relativistic electrons and magnetic fields, is equal to the total pressure (full line), thermal 
pressure (dashed-dotted line) and the kinetic pressure (dotted line). 

We observe that in many places the contribution of the kinetic pressure can be very 
important. This is a feature that is not being considered by the majority of works in the 
literature. It is very interesting to note that the surface brightness profile indicated by the 
total pressure (kinetic plus thermal) is very similar to that of the 151MHz data (see Figure 
11). The surface brightness profile indicated by the thermal pressure alone is very similar 
to that of the 1345 MHz data (see Figure 11). This suggests that different forms of energy 
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within the bridge may enhance radio emission at different frequencies differently. It is easy 
to see how this could occur since the radio emission at different frequencies is produced by 
relativistic electrons with different Lorentz factors. Thus, it would only require that the 
thermal pressure produces an underlying relativistic electron population with a power law 
distribution at all Lorentz factors, while the kinetic pressure boosts additional relativistic 
electrons to Lorentz factors that produce the 151 MHz emission but not the 1345 MHz 
emission. For example, through hydromagnetic processes the kinetic pressure could boost 
relativistic electrons with lower Lorentz factor up to the Lorentz factors needed to produce 
151 MHz emission. This could help to explain some of the changes in the radio spectrum 
seen across the source. 

To investigate this further, we compare the average thermal pressure, and the total 
average pressure (thermal plus kinetic pressures) obtained from the numerical simulation 
with the average pressure obtained from the 151 MHz and 1345 MHz data (Figure 25). To 
do this, the physical value of the units on the pressure obtained from the simulation were 
investigated, and were determined to be 1.596 x lO^^na dyn/cm , where Ua is the ambient 
gas density. A value of Ua = 10~^ cm~^ was adopted since this is consistent with X-ray 
measurements of the gas in the vicinity of the radio source (e.g.. Smith et al. 2002). The 
behavior of the pressures in the simulations roughly tracks that seen in the data. Of course, 
it is possible that the value of b in the hot spot region is closer to unity (e.g., Hardcastle 
et al. 2004; Donahue, Daly, & Horner 2003; Hardcastle et al. 2002; Wilson, Young, & 
Shopbell 2000), while that in the bridge is ~ 0.25. (e.g., Carilli et al. 1991; Perley & 
Taylor 1991; Wellman et al. 1997b). The value of b does not affect the pressures determined 
using the numerical code (which only depend on the scaling factor Ha), but it does affect 
the pressure determined from the data. If b in the hot spot in closer to unity, then the 
pressure determined using the data will be lower in this region that the values shown. This 
would improve the fit to the 1345 MHz data, and worsen the match with the 151 MHz data. 
Another effect to consider is the efficiency factor (converting gas kinetic into relativistic 
energy) used in converting surface brightness into pressure, which could vary with location 
within the source or could have an energy dependence, resulting in a frequency dependence. 
In addition, consider particle re-acceleration within the bridge region or a backflow could be 
important in the source. In view of these factors, the rough overall agreement between the 
data and the simulation is encouraging. 

Finally, in Fig. 26 we compare the thermal pressure determined from the simulation at 
several different times with that estimated from the 151 MHz empirical data for the two sides 
of Cygnus A, each normalized to match at the peak of the hot spot. Sources are dynamic, 
constantly changing with time, so we consider snapshot corresponding to panels c, d, and e of 
Figure 24 re-scaled to have a length of 60 kpc. The dynamic nature of the source is obvious 
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from these snapshots. The agreement is better than for the surface brightness distribution 
and this is due to the fact that the simulated pressure profile does not depend upon the 
bridge width. As we shall see below, the simulated and observed width do not agree equally 
well over the entire bridge extent. We see that the pressure profile varies considerably with 
time over the course of the simulation. This implies that there is no "steady-state" pressure 
distribution in the lobes. 



3.4. Bridge Width 

We compare the "true" width of the bridge in the simulation (Figure 20a) with the 
FWHM and second moment of the synthetic surface brightness distribution (Figure 27), 
where only the thermal component of the surface brightness was considered. We find that 
while the second moment gives a good fit to the true cocoon width, the Gaussian FWHM 
deviates significantly. (Note: this is not seen in the data, where the second moment and 
FWHM of the Gaussian fits track each other quite well Figure 10). The main departure in 
the simulation occurs between 20 and 30 kpc from the hot spot. This can be attributed to 
the fact that the cross-sectional slices of the surface brightness distribution in this region 
exhibits a double peak as can be seen in Fig. 23. This double peak implies that the gas 
pressure distribution inside the simulated cocoon is higher near its edges, at least in the 
region in question (20 — 30 kpc). 

The ratio ( = W0/W2 between the two quantities is in the range ~ 1.2 — 1.8 with an 
average value 1.39 ±0.07. For a perfect Gaussian distribution, the ratio ( = V21n2 ~ 1.177. 
To compare this with the empirical data shown in Fig. 10 we note that in this figure we have 
plotted 2^ = W^G'/(second moment). Since the value of the empirical data is in the range 
2 — 3 we conclude that the results of the simulations 2( ~ 2.4 — 3.6 are in fair agreement 
with the data. 

Figure 4 shows a comparison between the observed FWHM of the radio source with the 
"true" width of the simulated source. We see that the simulation shows the same general 
behavior as the radio source. 



3.4-1 ■ Emissivity Distribution 

As mentioned above, the value of ratio ( = Wg/W2 will depend on the emissivity 
distribution across the cocoon. In this sense, the ratio may be used as a diagnostic of the 
distribution of emissivity. 
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We have thus investigated this by calculating the ratio ( for a simple model cocoon 
whose radius varies smoothly according to 

where D is the distance from the hot spot and Re is in kpc. We studied five cases with 
different pressure profiles P{r). The radial dependence of P is shown in Table 1. In all cases 
P = for r > Re. 

We calculated the surface brightness using the same beam size as the empirical data 
and calculated the second moment and fit a Gaussian to the surface brightness profile per- 
pendicular to the cylinder axis in order to determine the ratio (. The values of W, Wq and 
W2 for four of the model bridges in Table 1 are shown in Fig. 29. We can clearly see how the 
relation between the true width W and the calculated ones Wq and W2 depends strongly on 
the pressure distribution across the cocoon. 

In Table 1 we show the average values of the Gaussian width, the second moment and 
(^ for the five distributions. We see that for Model A both the Gaussian fit and the second 
moment are about 55% to 70% of the true width, and thus need to be normalized to be 
used as a direct measure of the bridge width. For the other models the Gaussian width 
increasingly departs from the true width as the peak of pressure distribution moves away 
from the source axis. On the other hand, the second moment gives a better estimate when 
the pressure distribution concentrates near the border of the cocoon. 

We notice that the range of values of the ratio C, of the Gaussian FWHM and the 
second moment width W2 for the five distributions is in agreement with the value from the 
numerical simulations (1.2 — 1.8). The average value observed in Cygnus A (~ 2.5 = 2 x 1.25) 
is compatible with an emissivity distribution concentrated near the source axis (Model A). In 
a few cases when it reaches higher values (~ 3 = 2 x 1.5), the pressure distribution in these 
regions must have two separate peaks as in Model E. In fact, if we look at the deconvolved 
emissivity distribution of Cygnus A at 151 MHz (Figs. 15, 16, 17, 18) we see that this seems 
indeed to be the case in several sections of the bridge. Thus, the ratio between the Gaussian 
width and the second moment can provide a diagnostic of the emissivity profile in the lobe. 
It may be necessary to run 3-D numerical simulations to properly "calibrate" this diagnostic 
ratio. 

We have seen that, for the numerical simulation, the second moment of the surface 
brightness distribution gives a more stable estimate of the bridge width in our model cocoon 
than the Gaussian fit (e.g.. Fig. 27). If this applies to real sources it can be used as a 
measure of its real width, provided we have a proper normalization factor. Our simulation 
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gives an average value for the ratio between the cocoon and Gaussian width j3g = W/Wq of 
I3g = 0.91 ±0.12; while the average ratio of the true width to the second moment (32 = W/W2 
is p2 = 1-26 ± 0.04. 

We note that (using ( = Wg/W2 as a diagnostic, Table 1) the range of values of ( found 
in Cygnus A indicates that, at least in some regions, the bridge contains high pressure edges. 
Carilli et al. (1991) suggest that an edge brightening of the lobes could account for the 
increase in the magnetic field from the center to the edge of the lobes needed to explain the 
estimated age gradient. If this is the case, the Gaussian fit overestimates the bridge width 
as in Models D and E (Table 1) and, again, the second moment should be used to calculate 
the bridge width. 



4. SUMMARY 

We present a comparison of radio observations of the archetypal powerful radio galaxy 
Cygnus A and 2-D numerical hydrodynamical simulations. We characterize some global 
trends in the observed radio properties and compare them with the properties of a simulated 
radio source. We have analyzed detailed radio images at 151 MHz (Leahy et al. 1989) 
and 1345 MHz (Carilli et al. 1991). We have determined the surface brightness, emissivity, 
Gaussian FWHM, first and second moment of cross-sectional slices of the source as a function 
of distance along the bridges. We have also estimated the average pressure and minimum 
energy magnetic field in the radio plasma as a function of distance along the bridges. We 
find the following trends in global properties. 

• The first moment of the brightness distribution perpendicular to the source axis is 
a measure of the surface brightness weighted center of the source. We find that the 
first moment wanders with a peak-to-peak amplitude of ~ 10 arcsec. The sense of the 
wander has the opposite sign on the two sides of the source, but there is not an exact 
correspondence. 

• The width of the surface brightness distribution perpendicular to the source axis can 
be well characterized by a Gaussian over most of the length of the source. The ratio of 
the Gaussian FWHM to the second moment is fairly constant along the source with an 
average value of about 2.5 indicating that they give roughly consistent measurements 
of the source width. The radio lobe increases in width with distance from the hot spot. 
Superimposed on this general expansion are several "wiggles." 

• The average surface brightness, estimated pressure, and estimated minimum energy B 
field decrease with distance from the hot spots. 
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• We find evidence for significant structure in the estimated cross-sectional slices of 
emissivity. 

Following the methods of Carvalho & O'Dea (2002a,b) we present the results of 2-D 
numerical hydrodynamical simulations of light (density contrast ( = 10~^), supersonic (jet 
Mach number M = 10) jets propagating in a constant density environment. The results 
from the simulations are the following. 

• Jets propagating in a constant density atmosphere will decelerate with time (see also 
Carvalho & O'Dea 2002a). Thus, the estimated dynamical age of the source will 
be greater than the actual age of the source. For a source similar to Cygnus A the 
difference will be about a factor of 2. 

• We have run the simulations with several different grid widths. In the largest grid, 
the entire source is contained on the grid, while in the smaller grids, the bow shock 
eventually leaves the grid laterally. We find that the details of the cocoon properties 
differ in these simulations, though the overall properties of the cocoons are similar. 

• The pressure profile varies considerably with time over the course of the simulation. 
This implies that there is no "steady-state" pressure distribution in the lobes. There 
is considerable time dependence in the variation of the average surface brightness in 
the lobe near the hot spots, possibly due to vortex shedding. 

• The kinetic pressure in the cocoon (due to large scale turbulence) can make a significant 
contribution to the total pressure of the cocoon. We speculate that the thermal pressure 
produces an underlying relativistic electron population with a power law distribution 
at all Lorentz factors, while the kinetic pressure adds "extra" relativistic electrons 
with the Lorentz factors that produce the 151 MHz emission and not the 1345 MHz 
emission. For example, through hydromagnetic processes the kinetic pressure could 
boost relativistic electrons with lower Lorentz factor up to the Lorentz factors needed 
to produce 151 MHz emission. This could help to explain some of the changes in the 
radio spectrum seen across the source. 

• The second moment gives an accurate representation of the "true" width of the sim- 
ulated source. The Gaussian FWHM tends to be about 40% larger than the true 
width and can be systematically in error if the surface brightness exhibits multiple 
peaks. We suggest that the ratio of the Gaussian FWHM to the second moment may 
be a diagnostic of the emissivity profile in the lobes. However, "calibration" of this 
diagnostic probably requires 3-D simulations. Applying this diagnostic reveals that 
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bridge contains high pressure edges. This is consistent with the suggestion of Carilh et 
al. (1991) that an edge brightening of the lobes could account for the increase in the 
magnetic field from the center to the edge of the lobes needed to explain the estimated 
age gradient. 

We compare our simulations with the properties of the observed radio sources in two 
complementary ways. We note that the pressure is directly determined in the simulation but 
is estimated from the radio data, while the surface brightness is directly determined from 
the radio data but is estimated from the simulations. 

The simulations can qualitatively reproduce the overall observed morphology and the 
behavior of the cross-sections in surface brightness, the decline in surface brightness with dis- 
tance from the hot spots, and the width of the lobes. This suggests that the 2-D simulations 
give a reasonable representation of the properties of Cygnus A. 

This is especially encouraging given that the outflow axis of the jets in Cygnus A may 
wobble while those in the numerical simulation are held steady, and given that the ambient 
gas density in the simulation is assumed to be constant. The good agreement between the 
simulation and the data suggest that these effects are relatively small for a source like Cygnus 
A. 
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Table 1. Comparison between calculated and true width for a model bridge 

Pressure Gaussian fit Second moment 
distribution < Wg/IRc > < W2/2R^ > <( = Wg/W2 > 
P(r) 
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l-t 


0.673 


0.550 


1.223 


B 
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1.075 


0.786 


1.367 
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(t)"; 


1.211 


0.848 


1.428 


D 


(t) 


1.301 


0.889 


1.463 


E 


(t) 


1.416 


0.939 


1.507 
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Fig. 1. — Examples of cross-sectional slices used to define the surface brightness profiles 
shown in Fig. 2 for the 151 MHz data (left) and the 1345 MHz data (right). The distance 
of the slice from the hot spot is indicated at the bottom on the slice. The images have been 
rotated slightly for ease of display and the left side is the south-east lobe and the right side 
is the north-west lobe. 
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Fig. 2. — Cross-sectional surface brightness profiles (solid lines) and Gaussian fits (dashed 
lines) for the Left side of Cygnus A at 151 MHz (left panel) and 1345 MHz (right panel). 
The distance of the slice from the hot spot is indicated for each slice. 



-25- 






-20-10 10 20 -20-10 10 20 -20-10 10 20 -20-10 10 20 
arcsec arcsec arcsec arcsec 



"«, 


L=45.0 


Dh 


i^-'lB.G 


^Hi 


-51,0 


^HS 


-54.0 


: 


n 




A 








n, 


Dfis 


,=57.0 














r 


.lIlijjI 


..ulllljjj 




iImmImmImmIi 




iImmImmImmIi 




ImmImmImmIi 



-20-10 10 20-20-10 10 20-20-10 10 20-20-10 10 20 
arcsec arcsec arcsec arcsec 



Fig. 3. — Continuation of the cross-sectional surface brightness profiles and Gaussian fits for 
the Left side of Cygnus A at 151 MHz (left panel) and 1345 MHz (right panel). 
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Fig. 4. — Cross-sectional surface brightness profiles (solid lines) and Gaussian fits (dashed 
lines) for the Right side of Cygnus A at 151 MHz (left panel) and 1345 MHz (right panel). 
The distance of the slice from the hot spot is indicated for each slice. 
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Fig. 5. — Continuation of the cross-sectional surface brightness profiles and Gaussian fits for 
the Right side of Cygnus A at 151 MHz (left panel) and 1345 MHz (right panel). 
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Fig. 6. — The full width at half max of the best fit Gaussian as a function of slice distance 
from the hot spot at 151 MHz (left panel) and 1345 MHz (right panel). The upper panels 
show results from the Left side of the source, while the lower panels show results from the 
Right side. 
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Fig. 7.— Ratio of the FWHM obtained at 151 MHz to that obtained at 1345 MHz. A more 
narrowly focused view is shown on the right panel. 
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Fig. 8. — The first moment of the cross-sectional surface brightness profiles as a function of 
distance from the hot spot. The 151 MHz data is shown on the left panels, while the 1345 
MHz data is shown on the right panels; the Left side of the source is shown in the upper 
panels, and the Right side is shown in the lower panels. 
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Fig. 9. — The second moment of the cross-sectional surface brightness profiles as a function 
of distance from the hot spot for the 151 MHz data (left panel) and the 1345 MHz data 
(right panel). The upper panels show results from the Left side of the source and the lower 
panels show results from the Right side of the source. 
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Fig. 10. — The ratio of the Gaussian FWHM to the second moment as a function of distance 
from the hotspot for the 151 MHz data (left panel) and the 1345 MHz data (right panel). 
Measurements from the left side of the source are shown in the upper panels; the average 
values of these are 2.59 ± 0.07 and 2.45 ± 0.03 for the 151 and 1345 MHz data respectively. 
Measurements from the right side are shown in the lower panels; the average values are 
2.49 ± 0.07 and 2.53 ± 0.03 for the 151 and 1345 MHz data respectively. 
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Fig. 11. — The average surface brightness in a cross-sectional surface brightness shce as a 
function of distance from the hot spot for the 151 MHz data (left panel) and the 1345 MHz 
data (right panel) obtained assuming the diameter of the bridge at each distance from the 
hot spot can be approximated by the FWHM of the best fit Gaussian. Results from the Left 
side of the source are shown in the upper panels, while those from the Right are shown in 
the lower panels. 
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Fig. 12. — The average pressure in a cross-sectional surface brightness shce as a function 
of distance from the hot spot for the 151 MHz data (left) and the 1345 MHz data (right). 
The pressure for the 151 MHz data has units of 1.2 x 10^^^ erg cm~^(1.336~^'^ + 6^), where 
b parameterizes the offset from minimum energy conditions. For minimum energy con- 
ditions, 6=1, and the normalization factor or pressure unit for the 151 MHz data is 
2.8 X 10~^^ erg cm"^, while for a field strength that is about 0.25 the minimum energy value, 
the normalization is 1.3 x 10"^'' erg cm~^. The normalization factor for the 1345 MHz data is 
j^Q-io gj,g cixi-3(^x 33^-i-5_|_^2^ Pqj, niinimum energy conditions the pressure unit for the 1345 
MHz data is 2.4 x 10~^° erg cm~^, while for b = 0.25, the pressure unit is 1.1 x 10~^ erg cm~^. 
Measurements from the Left side of the source are shown in the upper panels, while those 
from the Right side are shown in the lower panels. 
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Fig. 13. — The average minimum energy magnetic field strength in a cross-sectional surface 
brightness slice as a function of distance from the hot spot for the 151 MHz data (left) and 
the 1345 MHz data (right). The units of the minimum magnetic field strength for the 151 
MHz data are 30 /iG, while those for the 1345 MHz data are 88 /iG. Estimates from the 
Left side of the source are shown in the upper panels, while those from the Right side of the 
source are shown in the lower panels. 
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Fig. 14. — Method of obtaining the emissivity as a function of radial distance from the shce 
center. 
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Fig. 15. — The radio emissivities (dashed hnes) and cross-sectional surface brightness profiles 
(solid lines) for the 151 MHz data (left panel) and the 1345 MHz data (right panel) for the 
Left side of Cygnus A. The units for the 151 MHz data are 1.8 x 10"^'* erg s~^ cm~^ Hz~^, 
while those for the 1345 MHz data are 9.0 x 10~^^ erg s~^ cm~'^ Hz~^. 
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Fig. 16. — Continuation of the radio emissivities and cross-sectional surface brightness pro- 
files for the Left side of Cygnus A. 
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Fig. 17. — The radio emissivities (dashed hnes) and cross-sectional surface brightness profiles 
(solid lines) for the 151 MHz data (left panel) and the 1345 MHz data (right panel) for the 
Right side of Cygnus A. The units for the 151 MHz data are 1.8 x 10"^'* erg s~^ cm~^ Hz~^, 
while those for the 1345 MHz data are 9.0 x 10~^^ erg s~^ cm~'^ Hz~^. 
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Fig. 18. — Continuation of the and radio emissivities and cross-sectional surface brightness 
profiles and radio emissivities for the Right side of Cygnus A. 
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Fig. 19. — Density contour map of a jet propagating in a constant density atmosphere for 
Run DOl (M = 10,7] = 2 X 10"^). The grid size is 66 x79 kpc. 
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Fig. 20. — (a) Contour of the bow shock, cocoon and jet for the simulation shown in Fig. 19 
(Run DOl). Panel (b) shows the longitudinal profile of the average cocoon pressure, cocoon 
lateral speed (c) and bow shock lateral speed (d). 
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Fig. 21. — Advance speed of the head (open circles) relative to the jet initial speed as a 
function of source size. The dashed line represents the ambient sound speed. 
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Fig. 22. — Synthetic radio image showing the time evolution of the numerical simulation of 
a jet (M = 10, ?7 = 2 X 10"*^) propagating in a constant density atmosphere. Each panel 
measures 144 x54 kpc and the synthetic data have been smoothed by a Gaussian with 
FWHM of 3 arcsec to match the 151 MHz data, as shown by the black circle in panel e. The 
age of the source corresponding to each panel is: (a) 1.35 Myr, (b) 4.53 Myr, (c) 6.17 Myr, 
(d) 8.89 Myr and (e) 9.15 Myr. One can clearly see the characteristic edge-brightened shape 
of powerful FR-II sources and an aspect ratio resembling that of Cygnus A. 
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Fig. 23. — Cross- sectional surface brightness profiles (solid lines) and Gaussian fits (dotted 
lines) for the numerical simulation (Run DOl). The distance between each panel is 3 kpc. 
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Fig. 24. — Average surface brightness at 151 MHz in a cross-sectional slice as a function 
of distance from the hot spot taking into account radiative losses. The radio emissivity is 
calculated supposing the relativistic electron density is proportional to the total pressure (full 
line), thermal pressure (dashed-dotted line) and the kinetic pressure (dotted line). Panels 
(a), (b) and (c) correspond to panels (c), (d) and (e) of Fig. 22, respectively. 
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Fig. 25. — Comparison between the logarithm of the observed and calculated average total 
(thermal plus kinetic - thick line) and thermal (thin line) pressures as a function of distance 
from the hot spot for 151 MHz (left panels) and 1345 MHz (right panels). The Left side of 
the source is shown in the upper panels, and the Right side is shown in the lower panels. The 
calculated pressure is for an ambient density Ua = 0.01 cm~^, and the units on the pressure 
are the same as in figure 12 for a value of the parameter b, defining the offset from minimum 
energy conditions, of 0.25. Thus, the unit on the 151 MHz pressure is 1.3 x 10~^° erg cm~^, 
while that on the 1345 MHz pressure is 1.1 x 10~^ erg cm~^. 
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Fig. 26. — Comparison between observed and calculated average thermal pressure as a 
function of distance from the hot spot, with the numerical result scaled to match the empirical 
result at the hot spot peak. The Left side of the source is shown in the upper panel, and the 
Right side is shown in the lower panel. The dashed line, dotted line and full line correspond 
to panels (c), (d) and (e) of Fig. 22, respectively. 
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Fig. 27. — Bridge width calculated from the Gaussian fit (FWHM) and the second moment 
compared with the full width of the cocoon measured from the simulations as a function of 
distance from the hotspot. Thick line is cocoon (W), thin line is the Gaussian fit (Wg) and 
dashed line is the second moment {W2)- 
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Fig. 28. — Comparison between the FWHM fit to the radio data with the "true" width of 
the simulated source as a function of distance from the hot spot for 151 MHz (left panels) 
and 1345 MHz (right panels). The Left side of the source is shown in the upper panels, and 
the Right side is shown in the lower panels. 
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Fig. 29. — Bridge width calculated from the Gaussian fit (FWHM) and the second moment 
compared with the full width of the model cocoon as a function of distance from the hotspot. 
Thick line is cocoon (W), thin line is the Gaussian fit (Wg) and dashed line is the second 
moment {W2)- The four panels represents the models in Table 1. 
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